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ABSTRACT 

Evidence for a second main-sequence turn-ofF in a deep cofour-magnitude diagram 
of NGC 1868 is presented. The data were obtained with HST/WFPC2 and reach 
down to 771555 — 25. Besides the usual r ~ 0.8 Gyr turn-ofF found in previous anal- 
yses, another possible turn-off is seen at 771555 — 21 {My — 2.5), which is consistent 
with an age of r ~ 3 Gyrs. This CMD feature stands out clearly especially when 
contaminating field LMC stars are statistically removed. The background subtracted 
CMD also visibly displays a red subgiant branch extending about 1.5 mag below the 
younger turn-ofF and the clump of red giants. The significance of the secondary turn-off 
in NGC 1868 was confirmed with Monte-Carlo simulations and bootstrapping tech- 
niques. Star-counts in selected regions in the cluster CMD indicate a mass ratio of old 

population/young population in the range 5% ^ Moid/Myoung ^ 12%, depending on 
the mass function slope. The existence of such a subpopulation in NGC 1868 is signif- 
icant even in the presence of uncertainties in background subtraction. The possibility 
that the secondary turn-off is associated with the field star population was examined 
by searching for similar features in CMDs of field stars. Statistically significant ex- 
cesses of stars redwards of the main-sequence were found in all such fields in the range 
20 ~ 771555 ^ 22. These however are much broader features that do not resemble the 
main-sequence termination of a single population. We also discuss other alternative 
explanations for the feature at 771555 ~ 21, such as unresolved binarism, peculiar stars 
or CMD discontinuities associated with the Bohm-Vitense gap. 

Key vifords: Magellanic Clouds: clusters; clusters individual: NGC 1868; stars: statis- 
tics. 



1 INTRODUCTION 

NGC 1868 is a rich Large Magellanic Cloud (LMC) clus- 
ter located at a = 5'' 14"" and 5 = -63° 57' (J2000), ap- 
proximately 6° away from the LMC's centre. Its age has 
been estimated both from ground-based photometry and 
spectroscopy, often yielding discrepant results covering the 
3 10* yrs ~ r ~ 10^ yrs range (Flower et al. 1980, Hodge 
1983, Bica & AUoin 1986, Chiosi et al. 1986, Olszewski et 
al. 1991, Corsi et al. 1994 and references therein). Most age 
estimates come from the interpretation of optical colour- 
magnitude diagrams (CMDs), usually based on the posi- 
tions of either the main-sequence turn-ofF (MSTO) or the 
clump of red giants (RC) or both. However, none of the 
CMDs available until now have been deep enough to al- 
low probing the main sequence at magnitudes as faint as 
V ~ 22 with small photometric errors (see the web page 
on www.ast.cam.ac.uk/STELLARPOPS/LMCdatabase for 



a detailed list of references on NGC 1868 and other rich 
LMC clusters in the HST Cycle 7 program GO7307). 

In this paper we investigate a possible second, fainter 
and therefore older, MSTO at V ~ 21, based on a deep 
CMD of NGC 1868 built from HST/WFPC2 data. This 
second population of stars may be the result of a strong 
interaction with another cluster, or of the capture of a lower- 
mass, older cluster by NGC 1868. In fact, there is evidence 
for a growing number of clusters in the LMC which may 
have undergone strong interactions or mergers (Kontizas et 
al. 1993, Dieball & Grebel 1998, Leon et al. 1999, Dieball, 
Grebel & Theis 2000). Confirmation of a merger event in 
a cluster's history, however, requires detailed photometric 
and/or spectroscopic data. Sagar et al (1991), based on a 
CMD of NGC 2214 with two apparent supergiant branches, 
suggested that it was made up of two distinct populations. 
However, Lee (1992), Bhatia & Piotto (1994) and Banks et 
al (1995), based on larger and more accurate photometry. 
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found a single population in the CMD of NGC 2214. Given 
that signatures of mergers or strong interactions remain ob- 
servable in clusters for at least 1 Gyr (de Oliveira, Bica & 
Dottori 2000), NGC 1868, if confirmed as a merger product, 
may be a good laboratory for studying the dynamical effects 
and the final products of such events. Hence, it is essential 
to investigate closely the possibility that it is made up of 
two distinct populations. 

Other explanations for a feature similar to a MSTO 
exist. Unresolved binaries, for instance, arc usually brighter 
and redder than the main sequence of single stars. Therefore, 
an enhanced fraction of unresolved binaries could result in 
such a feature. Sudden changes in stellar structure, such as 
the onset of convcctivc envelopes, may lead to CMD features 
such as the Bohm-Vitense gap, that may also mimic a MSTO 
(Bohm-Vitense 1970). These possibilities are also considered 
in the analysis of our deep NGC 1868 CMD, as well as of 
other similar CMDs. 

The paper outline is as follows; in §2 we describe the 
HST/WFPC2 data used, both on NGC 1868 and on control 
field areas. We also discuss the process of removing contam- 
ination by field stars from the on-cluster CMD. In §3 we 
test the statistical significance of the candidate secondary 
MSTO by means of bootstrapping realizations. In §4 we 
use isochrones and star counts in selected CMD areas to 
extract useful information about this presumed secondary 
population in NGC 1868. In §5 we explore the possibility 
that the candidate MSTO is related to field stars, making 
use of CMDs obtained in control field areas. We also explore 
the alternative mechanisms that could yield to features sim- 
ilar to a MSTO at that position. Our main conclusions are 
presented in §6. 



2 THE DATA 

The WFPC2 on-cluster data are described in more detail 
in Santiago et al. (2001). The off-cluster data used for field 
star subtraction were reduced and analyzed by Castro et al. 
(2001). In both cases the images were combined and cali- 
brated using the standard pipeline procedure (Holtzman et 
al. 1995a,b). Photometry was carried out using DAOPHOT 
tasks and is described in detail by Santiago et al (2001) and 
Castro et al (2001). Zero-points, GTE and aperture correc- 
tions were also applied, again following standard procedures 
(Santiago et al 2001, de Grijs et al 2002a,b). 

Figure 1 shows the resulting CMDs. The left panel 
shows stars belonging to the on-cluster images. Two such 
images wore taken, one with the Planetary Camera (PC) at 
the centre of NGC 1868 and the other at its half-light radius. 
They arc the CEN and HALF fields as defined by Santiago 
et al. (2001). The CMD shown represents the final sample, 
including stars from both HALF and CEN fields and with 
no repeats. The total on-cluster solid angle is approximately 
7.5 arcmin. As explained in Santiago et al. (2001), the stars 
in the region common to the HALF and CEN fields had two 
independent photometric measurements and were thus used 
to determine uncertainties as a function of 771555 magnitudes 
and m555 — nisu colour. A total of 731 stars were found in 
the overlap region between the HALF and CEN fields. These 
empirical error determinations are an essential part of our 
upcoming analysis. The uncertainties were computed in bins 



0.5 mag wide in both rnsss and nisu. In each magnitude bin 
we computed the mean value and standard deviation (a) of 
the distribution of magnitude differences. The typical un- 
certainty of a single measurement was then assumed to be 
a I \/2. The photometric uncertainty analysis is described in 
more detail in Kerber et al (2002). 

The right panel in Figure 1 shows stars in the off-cluster 
area, located 7.3' away from the centre of NGC 1868 and pre- 
viously studied by Castro et al (2001). This off-cluster field 
corresponds to a single WFPC2 field, covering about 5 ar- 
cmin. Both CMDs include only objects that were classified 
as stellar sources in the classification schemes presented in 
Santiago et al (2001) and Castro et al (2001). No additional 
cleaning of remaining spurious objects or background galax- 
ies was made, resulting in some objects located well away 
from the cluster main-sequence (MS), the red giant branch 
(RGB) and the RC. Most of them should be contaminating 
unresolved background galaxies or faint stars in our Galaxy. 

The on-cluster stars display a clear MS that terminates 
at m555 ~ 19.2. Previous works find the MS termination at 
V = 19.2-19.3 (Corsi et al 1994, Brocato et al 2001). This is 
also roughly the magnitude of the RC (771555 — 77J.814 — 0.8), 
which contains He burning stars and is largely dominated 
by the cluster stars. In both panels, a clear RGB stretches 
downwards from the RC position to the bottom of the sub- 
giant branch (SGB) at 771555 — 22, 771555 — mgi4 ~ 0.8. These 
stars mostly belong to the old (r ~ 10 Gyrs) LMC field 
star population (Holtzman et al. 1997, Castro et al. 2001, 
Smecker-Hane et al 2002). Besides the main cluster and old 
field turn-offs, there are two additional features that may 
be interpreted as MSTOs in the CMD: one has mostly field 
stars with 771555 ~ 20.3 and had already been identified by 
Castro et al. (2001) as a r ~ 2 Gyrs field population; the 
other is fainter, at 777555 — 21, 0.4 < 777555 — 771814 < 0.5 and 
consists entirely of stars in the on-cluster data. It is this lat- 
est feature that we concentrate on, as a candidate NGC 1868 
subpopulation. The MSTOs described here are indicated in 
Figure 1. 



2.1 Removing contaminating field LMC stars 

In order to clean the left panel of Figure 1 from background 
field contamination, we statistically remove the control off- 
cluster CMD from the on-cluster one. Several procedures for 
doing so have been considered. One method is to match the 
off-cluster stars to corresponding stars in the on-cluster im- 
age using the estimated probabilities that two stars could 
bo independent photometric measurements of each other. In 
order to estimate these matching probabilities, we use the 
empirically determined photometric uncertainties in both 
images and assume a Gaussian error distribution. Cluster 
stars are then randomly removed according to their prob- 
ability of matching any of the off-cluster stars. In this ap- 
proach a reliable estimate of the photometric uncertainties is 
very important. As pointed out previously, we used the inde- 
pendent magnitude measurements for the stars in common 
between the CEN and HALF fields to empirically estimate 
standard deviations (a) in the error distribution at different 
magnitude bins (Santiago et al 2001, Kerber et al 2002). 

An alternative way of removing contaminating field 
stars is to bin both on-cluster and off-cluster data in mag- 
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Figure 1. On-clustcr (left panel) and ofT-clustcr (right panel) colour-magnitude diagrams. The large trapezium box indieates the locus 
of evolved stars, whose field star subtraction was carried out separately from that of MS stars. Empirically determined photometric 
uncertainties are shown on the left panel. 



nitude (or in magnitude and colour) and to subtract the 
histogram of the latter from that of the former. At each 
bin, a subset of the on-cluster stars, corresponding to the 
subtracted histogram, is randomly selected. The two ap- 
proaches have been tested and yield similar results. As our 
main concern is to assess the existence of a second, older 
population superposed to the dominant NGC 1868 popula- 
tion, the branch of evolved stars is critical, since it provides 
the best opportunity to single out stars from each individual 
population. Therefore, field subtraction was carried out sep- 
arately for the MS and for the SGB/RGB locus. This latter 
is marked on Figure 1. 

A more detailed discussion about the issue of field stars 
removal is presented in an upcoming paper (Kerber et al. 
2002). 

The cleaned CMD is shown in Figure 2. The CMD is 
now cut-off at the faint end (msss > 24) as this region is not 
relevant to the current analysis. Below this limit, sampling 
incompleteness would only further complicate field star sub- 
traction. We also cut it at bright magnitudes (msss < 19) 
to avoid regions where saturation effects start to take place 
(m555 ^ 19 and msSS ~ 17.8 for the HALF and GEN fields, 
respectively). The statistical removal of field stars hcis de- 
pleted the subgiant branch in the range 21 ~ msss ~ 22. A 
residual number of SGB stars brighter than msGs — 21 but 
much fainter than the cluster RC remains. As this excess of 
SGB stars reltive to the field is located at brighter magni- 
tudes than the candidate second MSTO, it is a first evidence 



for the reality of this feature. As for the secondary MSTO 
itself, at 771555 — 21, it was left untouched, which reflects 

the absence of field stars close to it. This again supports the 
reahty of the second population in NGC 1868. 



3 STATISTICAL SIGNIFICANCE OF THE 
CANDIDATE SECONDARY MSTO 

In this section we will assess whether the feature we tenta- 
tively identify as a secondary MSTO in the GMD of NGG 
1868 is statistically significant. More specifically we address 
the question of whether a similar feature could originate 
from a single main sequence through a random realization 
of photometric errors in the data. We explore this possi- 
bility by means of Monte-Carlo realizations of the cleaned 
NGC 1868 CMD. In each realization we collapse the cluster 
MS onto a fiducial line and redistribute the data, using the 
measured photometric uncertainties and assuming Gaussian 
error distributions in magnitude and colour. 

The main sequence fiducial line is determined by tak- 
ing the median 771555 — ^77814 colour at different magnitude 
bins. We apply a Scr clipping to the colour distribution in or- 
der to eliminate outliers and iterate until convergence. We 
initially define this MS line in bins spaced by 0.2 mag. A 
cubic spline interpolation is then used to fill this line in a 
much narrower (A771-555 = 0.02) binning. Each true MS star 
is then randomly assigned to one of the points along the MS 
fiducial line, using a two-dimensional Gaussian probability 
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Figure 2. Colour-magnitude diagram resulting from subtra<;ting field stars from the on cluster data. 



function whose standard deviations are given by the star's 
measured magnitude and colour uncertainties. We then re- 
distribute the data back onto the CMD plane, the magni- 
tude and colour of each artificial star again being a random 
realization of the same Gaussian distribution of errors. 

One hundred such Monte-Carlo simulations were car- 
ried out. A subset with 11 of the resulting CMDs is shown 
in Figure 3. The actual data are shown in the lower-right 
panel for comparison. The MS fiducial line is shown as a 
white line for guidance. The artificial CMDs have similar 
width as the data but are unable to reproduce the appar- 
ent turn-off at ms,-,,-, ~ 21, or the residual of the old field 
population turn-off (771555 22), whose stars are farther to 
the red of the MS and much more clumped together than 
any set of points in any of the simulated CMDs. The sim- 
ulations also fail to reproduce the main cluster MSTO at 
777555 — 19 and 70.555 — 70,814 — 0.3. This proves that scat- 
ter due to photometric errors alone cannot account for these 
features seen in the data. Notice that these experiments do 
not explicitly account for the effect of unresolved binaries 
on the CMD. The fact that the distribution of true stars is 
skewed redwards from the MS fiducial line is likely caused 
by unresolved binaries. Incorporation of unresolved binaries 
would help spread out the simulated stars towards redder 
colours, but is unlikely to reproduce the features mentioned 
above, especially if we consider that part of the effect is al- 
ready incorporated in the position of the MS fiducial line. 
We discuss the effect of binaries in more detail in 85. 2. 



4 THE CANDIDATE NGC 1868 SECOND 
POPULATION 

Assuming that the second MSTO at 777555 ~ 21 is real and 
belongs to the cluster, one wishes to determine what are the 
main characteristics of this population in NGC 1868. 

In Figure 4 we again show the cleaned NGC 1868 CMD. 
Superimposed to the data we show three Padova isochrones 
(Girardi et al. 2000). In plotting the isochrones, we assume a 
distance modulus of m — M = 18.5 to the LMC (Panagia et 
al. 1991). The data were extinction corrected assuming E(B- 
V) =0.04 as indicated. Conversion from E(B-V) to E(555- 
814) was done as described by Holtzman et al (1995a,b). 
The [Fe/H] and age (in Myrs) of each isochrone are shown 
in the figure. The [-0.4,800] isochrone fits the main NGC 
1868 population. The other isochrones are attempts to fit 
the second turn-off and SGB. The younger and more metal- 
rich isochrone with [Fe/H] = -0.4 and an age of 2500 Myrs 
provides the best fit to both turn-off and SGB regions. 

Using star counts and this best fitting isochrones, we 
can estimate the mass fraction of this candidate older NGC 
1868 population relative to the dominant and younger one. 
There are 11 stars in the range 2.5 ^ M555 ~ 1.5 
(21.0 ~ 70555 ~ 20.0) and located along the SGB. Be- 
sides these, 13 stars are located within the CMD area whose 
limits are: 20.85 < 771555 < 21.3, 777555 — mgu < 0.6 and 
TiBBB < 5.333 (771555 — 777814) + 18.55. This corresponds to 
the region where the second MSTO is clearly detached from 
the MS. Assuming that all such Nou = 24 stars are in fact 
associated with an older population belonging to NGC 1868 
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Figure 3. A sample of artificial CMDs built from the NGC 1868 MS fiducial line (white line) and error distribution. The lower-right 
panel shows the actual NGC 1868 data, for comparison 
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Figure 4. Field subtracted cluster colour-magnitude diagram. Some Padova isochrones are superimposed to the data; their [Fe/H] and 
age (in Myrs) are indicated in the figure, as well as the adopted reddening value. 



and using the best fitting isochrone, we infer the masses 
that correspond to the hmiting range in absolute magni- 
tudes: they ~ 1.36 rriQ for the basis of MSTO 
position and rrimax — 1-49 rriQ for Afsr,,-, = 1.5 (ms,-,,-, = 20) 
along the SBG. As for the dominant cluster population we 



find N-u. 



= 732 stars along the cluster MS with Msss < 2 



(m555 < 20.5) and fainter than the MS termination (we cut 
it at m555 = 19). These limits correspond to 1.63 ~ m/m© 

^ 2.10 for the best fitting isochrone to the NGC 1868 upper 

MS ([Fc/H],agc(Myrs) = -0.4,800). Assuming the cluster's 
present day mass function (PDMF) to be a power-law with 
fixed slope a within the mass ranges considered, the mass 
ratio will be given by 



Mou 



No. 



''max,young 



rmn,youngi 



My 



7nin,oldi 



Assuming a SaJpeter value for the MF (a = 2.35), 
and the limiting mass values quoted above for each of the 



subpopulation, we then infer Moid/ My 



= 0.06. For 



1.5 wo have Moid/Myour. 



0.08 and for a 



3.2, 



Mold/ Myoung ~ 0.05. The mass ratios are also fairly in- 
sensitive to adopting the ([Fe/H],age(Myrs)) = (-0.7,3500) 
Padova isochrone for the secondary population; in this case, 
the old population mass range would be 1.22 < m/rriQ < 
1.29 and Moid/ Myoung would change by no more than 50% 
for any choice of a. The mass ratio estimates are summa- 
rized in Table 1. 

In addition to the PDMF slope and the model mass- 



Table 1. Mass ratios for the two subpopulations of NGC 1868 

stars for different choices of PDMF slope and assumptions for the 
age, metallicity and mass range of evolved stars in the secondary 
population 



isochrone 


SGB mass range (solar) 


a 


Mold/^yonng 


[-0.4,2500] 


1.36/1.49 


2.35 


0.06 


[-0.4, 2500] 


1.36/1.49 


1.50 


0.08 


[-0.4, 2500] 


1.36/1.49 


3.20 


0.05 


[-0.7, 3500] 


1.22/1.29 


2.35 


0.09 


[-0.7, 3500] 


1.22/1.29 


1.50 


0.12 


[-0.7, 3500] 


1.22/1.29 


3.20 


0.06 



luminosity relation, another source of uncertainty in the 
mass ratio quoted above is the residual contamination by 
field stars to the Noid counts in the SGB region. Previous 
experiments with the field star subtraction (sec §2.1) meth- 
ods have led to variations of up to 35% iu the resulting value 
of Noid- If we assume this to be the uncertainty caused by 
field contamination, a similar relative error is expected to 
propagate into Moid /Myoung- Notice that this is a conser- 
vative reasoning, since residual contamination by field LMC 
stars should also affect the younger population. We should 
also point out that this 35% uncertainty is larger than the 
Poisson fluctuation expected for the Noid counts. 

An important question is how the stars in this candi- 
date secondary population are distributed throughout the 
cluster. Differences in spatial distribution between the two 
populations could indicate a recent merger, in which the two 
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Figure 5. Dots: projected distribution of stars located at the CMD position of tlie secondary MSTO. The contours are circles centred 
on the cluster. The numbers shown vertically are the percentage of all cluster stars interior to each circle. 



systems have not yet attained dynamical equilibrium. Fig- 
ure 5 shows the projected distribution of the stars located 
at the second turn-off; they are the same 13 stars used in es- 
timating the Mold/ Myoung ratlo earlier in this section. The 
11 more evolved stars along the SBG were left out of Fig- 
ure 5 because they represent just a statistical excess relative 
to field counts; hence, they are not necessarily true mem- 
bers of the older population. The contours represent circles 
around the cluster centre and the attached labels show the 
percentage of all cluster stars located inside them. We note 
that, within the fluctuations expected by the small numbers, 
the projected distribution of these stars is consistent with 
that of the entire cluster. 



5 ALTERNATIVE INTERPRETATIONS 

We have so far been interpreting the CMD feature at msss ~ 
21 and msss — rrian ~ 0.4 as a second MSTO associated 
with NGC 1868. In this section we discuss other possible 
interpretations. 

5.1 Is it a population associated with the LMC 
field? 

One obvious alternative is that the feature in the NGC 1868 
CMD is in fact a MSTO, but associated with stars belonging 
to the general LMC field population rather than to the clus- 
ter itself. Notice however that this is unlikely, since most con- 
taminating field stars have been statistically removed from 



the on-cluster CMD. Furthermore, Figure 5 shows that the 
stars located at the position of the second turn-off are con- 
centrated towards the cluster centre, a strong evidence in 

favour of a cluster origin for them. On the other hand, the 
star formation history in the LMC is known to be complex, 
especially in the last few Gyrs, when the bar was formed 
(Gallagher et al 1996, Elson ct al 1997, Smeckcr-Hane et al 
2002). Thus, small scale variations in the field CMDs are 
possible, making field subtraction more uncertain. 

One way to investigate possible MSTOs in field popu- 
lations is to search for similar features in the 7 CMDs of 

field stars studied by Castro et al (2001). These are fields 
located 7.3' from the target clusters of the GO7307 HST 
project. In fact, Castro et al have visually identified possi- 
ble turn-offs, with ages in the range 2 ~ r ~ 4 Gyrs, in 
several of these fields. In addition to those, turn-offs associ- 
ated with an old (> 10 Gyrs) population were found in all 
field CMDs. In some cases, visual inspection of the CMDs 
revealed a broadening in the MS for msss ~ 21.5, indicative 
of continuous star formation in the LMC. 

In order to quantify these visual impressions we tried to 
identify features in the field CMDs containing stars that are 
highly detached from the MS line, similar to the presumed 
NGC 1868 MSTO. One obvious such locus is the SGB/RGB 
region, which is always present in the CMDs of field LMC 
stars. It thus has to be eliminated from the analysis a priori. 
We must also avoid contamination from background galaxies 
and faint stars belonging to the Galaxy. Unresolved galax- 
ies should be limited to faint magnitudes, msss ~ 22.5, 
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Figure 6. CMDs of field LMC stars with ttie 5% tiigtiest values of rigig as defined in the text. Each panel corresponds to one particular 
field in the LMC as indicated. The central solid line is the MS fiducial line and the dashed lines correspond to its ± 5a deviations. 
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Figure 7. Frequency distribution of the stars shown in the previous figure as a function of msss magnitude. Again each panel corresponds 
to one of the LMC fields, as indicated. The histogram shown as a dotted line corresponds to the frequency distribution of all MS stars 
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as brighter ones are visibly extended in our WFPC2 im- 
ages. From counts of faint compact galaxies, one expects 
~ 50 contaminating galaxies within the 22.5 < msss < 25 
range in our on-cluster field (Abraham et al 1996). As for 
foreground stars, a similar number is expected in the entire 
observed CMD (Santiago et al 1996). Together both types 
of contaminating sources contribute with a few percent of 
the total sample and will be spread out in the CMD. 

We therefore proceeded as follows. For each field a fidu- 
cial line representing the CMD MS was defined in the same 
way as described in §3. A low-order polynomial was fit to 
the MS fiducial line of most fields in order to smooth out the 
wiggles caused by noise in the median 771555 — »7i8i4 value at 
each m555 bin. For NGC 1805 and NGC 1818, which have 
much larger numbers of MS stars than the other fields, espe- 
cially at the bright end (20.5 < msss < 19), the raw MS line 
was used. Evolved stars, possible background galaxies and 
faint stars belonging to the Galaxy wore then eliminated by 
cutting out all CMD objects beyond ±5 a from the MS line, 
where a is the empirically determined standard deviation 
in the 771555 — mau colour distribution at each magnitude. 
For the remaining stars, we computed the number of stan- 
dard deviations by which each star is detached from the MS 
fiducial line, Usig = A(?7i555 — mgii) fa- 
in Figure 6, we show the stars whose risig values fall at 
the 95% position or beyond in each LMC field. The num- 
ber of stars in each panel varies from 30 to 130. These stars 
would bo > 2 cr (and < 5 a) events of a Gaussian error 
distribution in colour. The MS fiducial line and the ±5 a 
lines are also shown for guidance. Assuming that the pho- 
tometric uncertainties have been adequately measured over 
the entire 771555 range and that these stars just reflect the 
high tail of the error distribution, we would expect them to 
follow the distribution of MS staxs. Therefore, the shape of 
the distribution of such stars may reveal features that are 
not accounted for by errors alone. As an example, a clumped 
distribution of such stars at some magnitude range, as com- 
pared to the smooth distribution of MS stars, would indicate 
the existence of features such as a MSTO. 

There is a common pattern in the distribution of these 
5% reddest stars in the 7 LMC fields. In all cases, the ma- 
jority of them are within the range 20 < 777555 < 22. This 
can be also seen by the frequency distribution of these stars 
as a function of 777555, which is shown in Figure 7. In all 
fields, largo Ti^^g stars show a distinct peak in their magni- 
tude histogram when compared to the distribution of nor- 
mal MS stars. This peak is usually broad and covers the 
777555 range mentioned above. Examples are the fields close to 
NGC 1805, NGC 1818 and NGC 1868. NGC 1805 and NGC 
1818 arc close to star formation regions and their neigh- 
bouring field includes stars of difi'erent magnitudes (masses) 
leaving the MS towards the RGB. As for the field close to 
NGC 1868, a 2 Gyr MSTO is seen in its CMD (Castro et al 
2001), contributing with stars brighter than 711555 ~ 20.5. 

In some cases, on the other hand, high Usig field stars 
show a narrow and very tall peak. This is the case of the 
field stars in the neighbourhood of NGC 1831, which display 
a distinct peak at 777555 — 21; this is consistent with the 
position of the ~ 4 Gyrs turii-ofl' identified in this field by 
Castro et al (2001). Other examples of high peaks in the 
distribution of detached field stars include the fields close to 
Hodge 14 and NGC 2209. These peaks are at least 0.5 mag 



fainter than that of NGC 1831 and are close to the basis of 
the SGB of the old LMC field populations. Inspection of the 
stars contributing to these peaks in fact confirms that they 
are dominated by SGB stars which were not successfully 
removed using the 5 cr colour cut-off. 

Figures 8 and 9 show, respectively, the CMD and 777555 
histogram for the stars in the on-cluster, background cleaned 
NGC 1868 data whose Usig values are at the 95% position or 
beyond. These figures follow the same conventions as Fig- 
ures 6 and 7. A total of 201 stars contribute to Figures 8 
and 9. The larger number of stars allowed a magnitude bin- 
ning in Figure 9 narrower than in Figure 7. Notice that the 
on-clustcr histogram of high risig stars as a function of mag- 
nitude now shows a large and narrow peak at a bright magni- 
tude, 177555 — 19.5, which is close to the dominant NGC 1868 
MSTO. A second peak is seen in the range 20.5 < 777555 < 22, 
similar to those found in the field CMDs. This latter peak 
is well centered where the candidate secondary MSTO was 
found. Apart from the bright peak, however, the distribu- 
tion of high 7isig stars as a function of 171555 in the cluster 
data is not markedly different from those of field stars. 

Thus, we conclude that a field origin for the second turn- 
off in NGC 1868 may not be ruled out. However, the differ- 
ences in shape between this feature and those known to be 
associated with field LMC stars, plus the fact that the NGC 
1868 CMD had previously been removed of field stars, argue 
against this interpretation. Furthermore, the suspected sec- 
ondary turn-off in NGC 1868 is detached enough from the 
MS that the two reddest stars contributing to it actually are 
located beyond the -|-5 a line shown in Figure 8 (see Fig- 
ure 2), and therefore do not contribute to the histogram in 
Figure 9. 

5.2 Selective unresolved binarism and Am stars 

As mentioned in §3, there are other causes, besides a MS 
turn-off, that could lead to a set of stars being significantly 
detached from the MS towards redder colours. 

Unresolved binarism is one such possibility. The pres- 
ence of a secondary star, whose fiux is added to the primary, 
will shift the position of the system towards brighter mag- 
nitudes in a CMD, relative to the position of the primary. If 
both are MS stars, the secondary will be redder, thus also 
affecting the colour of the system. Unresolved binaries are 
certainly present in any CMD and should be distributed all 
over the MS and among evolved stars. Their main effect is to 
make the MS broader with a red tail in the colour distribu- 
tion. The feature seen in Figures 1 and 2, at 777555 — 21, could 
be explained only if there is a specific increase in the fraction 
of cluster binaries at an absolute magnitude M555 ~ 2.5. 

The apparent second MSTO occurs at the luminosity 
typical of peculiar A stars (Ap, Am, etc). Due to their chem- 
ical anomalies. Am stars in particular tend to be redder than 
normal A stars by about 0.05 mag as a result of increased 
lino blanketing. This effect alone is too small to account for 
the feature in question but certainly goes in the right direc- 
tion. As many Am stars are also known to be binaries (Car- 
quillat et al 2001, Debernardi et al 2001), a combination of 
blanketing and enhanced binarism might accommodate the 
observed feature in the CMD of NGC 1868. However, more 
stringent constraints on the frequency of Am staxs and on 
their binary fraction must be placed in order to test this 
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Figure 8. CMDs of the NGC 1868 stars with tlic 5% highest values of n^jg as defined in the text. The central solid line is the MS 
fiducial line and the dashed lines correspond to its ± 5cr deviations. 



possibility. Another problem with this interpretation is that 
global binary fraction estimates tend to be smaller in rich 
clusters than in the general field (Mayor et al 1996, Cote et 
al 1996, Elson et al 1998). 



5.3 The Bohm-Vitense gap 

The onset of a convcctive envelope is known to occur for 
MS stars of spectral type A or later. This change in en- 
ergy transport is expected to occur abruptly for stars with 
Te// ~ 7500 K, or (S-V) ~ 0.22 (Bohm-Vitense 1958). As 
a result, stars cooler than this critical effective temperature 
would become redder than the slightly hotter ones with fully 
radiative envelopes (Bohm-Vitense 1970). The resulting dis- 
continuity in the CMD is the so-called Bohm-Vitense gap 
and is thought to have an amplitude of A(_B — 1/) ~ 0.10. 
Recent surveys combining precise photometry and distance 
measurements from astrometry have shown that a gap does 
exist and is located not too fax from where expected by the- 
ory. But the evidence is often weak or subject to selection 
effects. Newberg & Yanny (1998) studied bright and nearby 
field stars with available photometry plus distances from 
Hipparcos and found a gap at 0.2 < {B — V) < 0.3 in the 
CMD of field stars of luminosity class V only. In other words, 
the gap is masked out by evolved stars contaminating the 
apparent MS in their CMD. Rachford & Caterna (2000) de- 
tected a gap at {B — V) ~ 0.35 in most nearby open clusters 
they studied, but they argue that this gap may be unrelated 
to the onset of envelope convection. 



We have carried out the same gap detection analysis as 
Rachford & Caterna (2000). We computed values as de- 
fined by those authors for gap candidates within the 0.15 < 
'"sr).-) — fnsiA < 0.60 colour range and with widths varying 
from 0.02 to 0.10 mag. Figure 10 shows the msss — m8i4 fre- 
quency distribution, not only for NGC 1868, but also for 3 
other rich LMC clusters in our sample, each one with more 
than 5000 stars in their CMD (Santiago et al 2001). The 
clusters are paired according to age, NGC 1805 and NGC 
1818 being younger (r < 10® yrs, Johnson et al 2001) than 
NGC 1831 and NGC 1868. The colour distributions for the 
older clusters is slightly flatter than those for the younger 
ones. Apparent gaps are seen at msss — mm — 0.40 (in all 
cases) and at msss — msu — 0.52 (except NGC 1831). A 
broader gap is also present at msss — tusia — 0.25 — 0.30 in 
NGC 1805 and NGC 1818 and possibly in NGC 1831 and 
NGC 1868. 

Figure 11 shows the gap center plotted against gap 
width for all our gap candidates with > 90% confidence 
(X^ > 2.8). Each panel corresponds to one of the 4 rich clus- 
ters discussed. The gaps at msss — m8i4 — 0.4 are significant 
in all clusters except NGC 1818. At msss — m8i4 — 0.52, only 
narrow and fairly unconspicuous gaps are confirmed. The 
bluer gap, at msss — m8i4 — 0.25, takes different widths and 
positions, being more significant in NGC 1805 and NGC 
1831. Hence, the main feature, clearly present in 3 rich clus- 
ters, is the gap at 0.35 < msss — m8i4 < 0.40. Notice 
that this colour coincides closely with that of the candi- 
date MSTO, adding support to the statistical reality of this 
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Figure 9. Solid line: Frequency distribution of the stars shown in the previous figure as a function of msss magnitude. The dotted Une 
corresponds to the frequency distribution of all MS stars in NGC 1868 



CMD feature. These gap limits correspond approximately to 
0.23 < B — V < 0.27, which is close to the expected position 
of the Bohm-Vitense gap and somewhat displaced from the 

gap position favoured by Rachford & Catcrna (2000). 

The presence of colour gaps in most of our rich LMC 
clusters calls for a more common mechanism than merger 
events. Our gap detection results, therefore, favour the pos- 
sibility that the Bohm-Vitense gap was detected for the first 
time in LMC clusters. On the other hand, our gap colour is 
inconsistent with the one found by de Bruijne et al (2000, 
2001) in their very accurate Hyades data. In fact, the very 
narrow CMDs by these authors provide the strongest lim- 
its on the gap position and associated colour discontinuity, 
which seems to be of A{B — V^) ~ 0.05. This is significantly 
smaller than the colour difference between MS stars and the 
secondary MSTO found in our NGC 1868 data. 



6 CONCLUSIONS 

In this paper we have shown possible evidence for a sec- 
ond main sequence turn-off in a deep WFPC2/HST colour- 
magnitude diagram of NGC 1868. This feature is clearly visi- 
ble in the cluster CMD, at msss ~ 21 and msss— m8i4 — 0.4, 
especially after contaminating field stars are removed. The 
presence of these stars, as well as of a residual SGB/RGB 
brighter than ms,-,,-, — 21, is evidence for a secondary cluster 
population. In fact, previous ground-based data also show a 
hint of this feature: the distribution of stars seen redwards 
of the main-sequence in the NGC 1868 CMD by Corsi et al. 



(1994), despite the large scatter, has an excess of stars at 
the same position. 

Assuming that the feature is really a MSTO associated 
with NGC 1868, isochrone fits yield an age of 2.5 — 3.5 Gyrs 
and a metallicity of —0.7 < [Fe/H] < —0.4 for this subpop- 
ulation. Using star counts in CMD regions where we expect 
only one cluster subpopulation to be present, we estimate 
the mass ratio of the older subpopulation relative to the 
younger to be in the range 0.05 ~ Moid/Myoung ~ 0.12. 
This range of values incorporates uncertainties in mass func- 
tion slope, in age and metallicity of the secondary popula- 
tion, as well as uncertainties in the statistical removal of 
contaminating field LMC stars. 

Even though the candidate second turn-off remains un- 
touched after field stars are statistically subtracted from the 
CMD of NGC 1868, the possibility that this feature is as- 
sociated with a field star population was investigated in de- 
tail and cannot be completely ruled out. CMDs of field stars 
from several different positions in the LMC show an excess of 
stars located in loci close to that of the NGC 1868 candidate 
turn-off. This excess is measured relative to the expected 
number of stars scattered to these loci due to photomet- 
ric errors. Most of the features in the field CMDs, however, 
are broad and visually less conspicuous than the apparent 
secondary turn-off in NGC 1868, thus being more consis- 
tent with periods of enhanced field star formation lasting 
for longer than 1 Gyr. Another argument against a field ori- 
gin for the candidate second turn-oflf is that its stars are 
concentrated towards the cluster centre. 
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Figure 10. Colour distributions for NGC 1805, NGC 1818 (upper panel), NGC 1831 and NGC 1868 (lower panel). 



Other possible explanations for the CMD feature have 
been explored. They include CMD spread due to unresolved 
binaries, which could mimic a MS turn-off if binarism is 
enhanced within a narrow range of primary star masses. It 
is interesting to notice that the mass range in question is 
close to that of Am stars, which are redder than normal A 
stars and for which a larger than usual binary fraction may 
exist. Estimates of binary fraction within globular clusters, 
however, yield smaller values than in the field. Besides, the 
binary fraction among Am stars is not yet well constrained. 
It is thus unlikely that the CMD feature discussed here is 
due to Am stars or unresolved binaries. 

A final and exciting alternative would be that the 

Bohm-Vitense gap was for the first time detected in a stel- 
lar population outside the Galaxy. The fact that gaps were 
found in the colour distribution of two other LMC clus- 
ters besides NGC 1868, all of them at m^bz — msu ^ 0.4 
{{B — V) ~ 0.25), calls for a more common mechanism than 
mergers. In fact, a merger of the type suggested by the pres- 
ence of a second and much fainter MSTO as discussed in this 
paper should be rare (Vallenari et al 1998, Dieball & Grebel 
2000). Most model predictions favour merging units that are 
nearly coeval and usually formed through encounters within 
the same giant molecular cloud, cither before or after they 
are fully formed clusters (Fujimoto & Kumai 1997, Efremov 
& Elmegreen 1998, Leon et al 1999). 

On the other hand, it is unclear to what extent the 
Bohm-Vitense gap may mimic a main sequence turn-off. The 
recent and precise CMD data on the Hyades by de Bruijne 
et al (2000, 2001) do not reveal strong turn-off like features 



associated with the Bohm-Vitense gap candidates found by 
those authors; the colour discontinuity in their CMD is of 
smaller amplitude than necessary to account for the NGC 
1868 feature studied here. Besides, the gap associated with 
our NGC 1868 CMD feature is bluer than the gaps found 
by de Bruijne et al. 

We should point out that the nature of the stars that 
make up the candidate MSTO in NGC 1868, which currently 
is uncertain and accountable for by at least two astrophys- 
ically interesting alternatives (merging and Bohm-Vitense 
gap), may be established by spectroscopic classification, a 
task likely to demand large, 8m class, telescopes. 
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